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ABSTRACT 

We study the star formation rate (SFR) as a function of environment for UV-selected Lyman break 
galaxies (LBGs) at redshift 3. From deep {^J-I,AB{sky) ~ 27.6) UBVI KPNO 4-m/MOSAIC images, 
covering a total of 0.90 deg^, we select 334 LBGs in slices I00/i“^ Mpc (co-moving) deep spanning the 
redshift range 2.9 < z < 3.4 based on Bayesian photometric redshifts that include the I magnitude as 
a prior. The slice width (100/i“^ Mpc) corresponds to the photometric redshift accuracy (A^ ^ 0.15). 
We used mock catalogs from the GIF2 cosmological simulations to show that this redshift resolution is 
sufficient to statistically differentiate the high-density regions from the low-density regions using E5, the 
projected density to the fifth nearest neighbor. These mock catalogs have a redshift depth of 110/i“^ Mpc, 
similar to our slice width. The large area of the MOSAIC images, ~ 40 x 40 Mpc (co-moving) per field, 
allows us to measure the SFR from the dust-corrected UV continuum as a function of E5. In contrast 
to low-redshift galaxies, we find that the SFR (or UV luminosity) of LBGs at z = 3 shows no detectable 
dependence on environment over 2 orders of magnitude in density. To test the significance of our result, 
we use Monte Carlo simulations (from the mock catalogs) and the same projected density estimators 
that we applied to our data. We find that we can reject the steep z = 0 SFR-density relation at the 
5—cr level. We conclude that the SFR-density relation at z = 3 must be at least 3.6 times flatter than 
it is locally, i.e. the SFR of LBGs is significantly less dependent on environment than the SFR of local 
star-forming galaxies. We find that the rest-frame UV colors are also independent of environment. 

Subject headings: cosmology: observations — galaxies: high-redshift — galaxies: evolution 


1. INTRODUCTION 

There is a long history of results showing that galaxy 
properties vary as a function of environments such as the 
morphology-density relation in clusters (Melnick & Sar¬ 
gent 1977; Dressier 1980; Postman & Geller 1984; Dressier 
et al. 1997; Treu et al. 2003) and its associated star forma¬ 
tion rate-density relation (e.g., Lewis et al. 2002). Lewis et 
al. (2002) used the Two-Degree Field Galaxy Redshift Sur¬ 
vey (2dFGRS; Colless et al. 2001) to show that the anti¬ 
correlation between star formation rate (SFR) and local 
projected density holds in clusters up to 2 virial radii from 
the cluster center. Not only is the mean luminosity a func¬ 
tion of environment (Blanton et al. 2003; Hogg et al. 2003; 
Balogh et al. 2004), but the luminosity function (LF) is 
also a continuous function of the local density (e.g. Cro¬ 
ton et al. 2004). 

The SFR also depends on the local density in groups 
and in the field. Gomez et al. (2003) showed similar re¬ 
sults using the Early Data Release (Stoughton & et ah, 
2002) of the Sloan Digital Sky Survey (SDSS), and they 
also showed that the SFR-density relation holds for field 
galaxies. Balogh et al. (2004) used both the 2dFGRS and 
the SDSS to show that the SFR (as measured by the Ha 
equivalent width) of field galaxies is strongly dependent on 
the local projected density. All these results point to the 
existence of physical mechanisms that quench star forma¬ 


tion as the local density increases, from the field to groups 
to clusters. 

At z = 0.4, the quenching of star-formation is appar¬ 
ently already in place. Using the wide-held (~ 30') im¬ 
ager SuprimeCam at the Subaru 8-m telescope, Kodama 
et al. (2001) selected members of the z = 0.41 cluster 
A851 using photometric redshifts and found that the col¬ 
ors of faint galaxies are strongly dependent on the lo¬ 
cal galaxy density. Follow-up narrow band Ha imaging 
of this cluster showed that the fraction of star-forming 
galaxies is a strong function of the local projected density 
far Mpc) from the cluster center (Kodama et al. 

2004). 

No observations of the SFR-density relation have yet 
been made at redshifts greater than one. Here, we used 
ground-based data from the Kitt Peak 4-m telescope htted 
with the MOSAIC wide-held imager to investigate whether 
or not the SFR-density relation already existed at redshift 
z = 3. We use photometric redshifts that are accurate to 
Sz ~ 0.15 or ~ 100 h~^ Mpc, comparable in co-moving 
distance to those of the COMBO-17 survey. 

We present our data and briehy discuss our selection of 
Lyman break galaxies (LBGs) in section 2. We then use 
mock catalogs from the G1F2 collaboration (Gao et al. 
2004) to show in section 3 that our photometric redshift 
accuracy is sufficient to statistically distinguish between 
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low and high density regions. We estimated the SFR as 
in section 4. Finally, our results are presented in section 5 
and conclusions in section 6. 

Throughout this Letter, we adopt Om = 0.3, Oa = 0.7 
and Hq = 100 h km s“^ Mpc“^. At redshift z ~ 3, 
H{z) ~ 4.46i7o, so 5z = 0.1 corresponds to 67 h~^ Mpc in 
co-moving coordinates. 

2. DATA 

The observations used here are described in detail in 
Bouche (2003) and Bouche & Lowenthal (2004, hereafter 
BL04). Briefly, we used the wide-held MOSAIC camera 
(Jacoby et al. 1998) at the Kitt Peak National Observatory 
4m telescope to image three helds in UBVI to a limiting 
magnitude ni,AB{sky) ~ 27.6. These helds were selected 
for the presence of a damped Lya absorber (DLA) at z ~ 3 
in our original DLA survey (Bouche & Lowenthal 2003, 
BL04). Given the camera’s held of view of 36'x36'(or 
~ 40 h~^ Mpc on a side), our survey covers a total of 0.90 
deg^. We selected ~3,000 LBGs within the redshift range 
2.8 < Zphot < 3.5. From Monte Garlo simulations, our 50% 
completeness level is Iab = 24.8 {TZ ~ 25), corresponding 
to ~ 0.7L*, using L*(z = 3) from Steidel et al. (1999). 

Photometric redshifts were obtained from the algorithm 
Hyperz (Bolzonella et al. 2000) coupled with the prior like¬ 
lihood distribution given by the magnitude mi of each 
galaxy (as in Benitez 2000). See Bouche (2003) and BL04 
for a detailed description of the tests we performed on the 
algorithm using spectroscopic redshifts in the Hubble Deep 
Field-North. Up to z < 6, the overall rms is 0.1-0.15, 
corresponding to ~ 67-100 h~^ Mpc, and A2/(l -|- Zgpec) 
is 0.06. 

In order to selected LBGs in a narrow redshift slice, we 
used the criterion: 

P(zo ± W,/2) > 0.5, (1) 

where zq is the slice redshift, 114 is the redshift slice width, 
and P{z) is the redshift probability distribution (BL04). 
We choose a redshift width of = 0.15 because, as dis¬ 
cussed in Bouche & Lowenthal (2003), it produces the 
largest sample in the smallest redshift slice, given the rms 
of the photo-z’s. At zq = 3, there are 68, 162, and 104 
galaxies in our three fields for a total of 334. 

3. CAN WE USE PHOTOMETRIC REDSHIFTS TO MEASURE 

Ss? 

Before investigating the galaxy properties as a function 
of environment, we used the GIF2 numerical simulations 
of galaxy formation (described in Gao et al. 2004 and De 
Lucia, Kauffmann, & White 2004) to investigate whether 
or not, with such a coarse redshift resolution, one can 
make reliable density estimations based on projected sur¬ 
face densities. These A-body simulations were run with 
400^ particles (each 1.73 x IO^Mq) over a periodic cube of 
110 h~^ Mpc on a side in a A dominated universe. From 
the output at z = 3 containing ~ 29,000 galaxies, we 
simulated an LBG sample by selecting 975 objects with 
mi < 25 and {B — V) < 0.2, where mi is the observed / 
magnitude given the absolute magnitude M/, the redshift 
of the simulation z = 3, and the AT-correction according 
to the type of spectral energy distribution of the galaxy. 


Fig. 1 shows the two-dimensional distribution (projected 
along the z-axis) of simulated LBGs (filled triangles) along 
with all the galaxies (points) in the simulated GIF2 vol¬ 
ume. One sees that the simulated LBGs do trace the ma¬ 
jor over-dense regions and filaments. Importantly, this 
‘slice’ width corresponds to our redshift resolution, i.e. 
~ 100 h-^ Mpc. 

Fig. 2 shows, for the GIF2 simulation, the two- 
dimensional surface density, E5, versus the three- 
dimensional true density ps. Both E5 and ps are com¬ 
puted from the distance to the fifth nearest neighbor. The 
horizontal error bars reflect the Icr rms to the mean. One 
order of magnitude change in corresponds to about 1.5 
orders of magnitude in p„. Thus, we conclude that the 
projected density gives a statistical measure of the true 
three-dimensional over-density even when it is projected 
over 110 h~^ Mpc. 

4. ESTIMATING STAR FORMATION RATES 

As reviewed by Kennicutt (1998), one can estimate the 
SFR from the rest-frame UV luminosity density Li, in the 
range 1500-2500 A: 

SFR(Mq yr“^) = 1.4 x 10“^®L^(ergs s“^ Hz“^) (2) 

for a Salpeter initial mass function, covering the range 0.1- 
100 Mq. This relation applies only to galaxies with con¬ 
tinuous star formation over time scales of 10® yr or longer. 
Pettini et al. (2001) have shown that for a small sample 
of LBGs, estimates of SFR from Eq. 2 are consistent with 
the SFR inferred from H/3 line fluxes. 

We estimated the rest-frame luminosity densities (in 
ergs s“^ Hz“^) of the 334 LBGs in our MOSAIG images 
from the standard expression = ioO'4 ("iab-i-48.6)^^2^ _I_ 
•Zphot) and from the luminosity distance. The luminosity 
was corrected for dust extinction. The amount of dust 
extinction {Ay) is estimated simultaneously with the es¬ 
timation of the photometric redshift Zphot- However, we 
found, using mock catalogs, that Ay recovered by Hyperz 
is biased: Ai/(true) = 0.51 Ay + 0.32. These mock cat¬ 
alogs were created from the same galaxy templates, with 
the same extinction law used in the redshift fitting and 
restricted to our redshift range of interest 2.85 < z < 3.25. 

We used the /-band (Ac = 8214 A) flux to estimate the 
SFR; this filter corresponds to a rest-frame wavelength of 
2000 A for z = 3 galaxies. We used the Calzetti et al. 
(2000) extinction curve ^ to estimate A °, and the dust 

corrected SFR (SFRg) follows from and 

Eq. 2. 

5. RESULTS AND DISCUSSIONS 

Eor galaxies at Zphot = 3 ± 0.075 selected with Eq. 1, 
Fig. 3 shows the dust-corrected star formation rate (SFRq) 
as a function of the normalized local density E5/ < S5 > 
where E5 is the surface density within the 5th nearest 
neighbor, i.e. E 5 = 6/(7rr|), and < S5 > is its median in 
each of our three MOSAIC fields. The local over-density 
(55 is E5/ < E5 > — 1. We use the normalized density in 
order to combine our several fields given that they do not 
contain the same mean number of LBGs per unit area. 
Error bars are shown for 20 randomly selected galaxies. 


This is the same extinction law that was used by Hyperz in determining the photometric redshifts. 
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There is no detectable difference between the distribution 
of SFRs of LBGs in low-density environments and those 
in high-density environments. 

We used the Kolmogorov-Smirnov (KS) test to deter¬ 
mine whether or not the galaxies in high density regions 
with S 5 > X are a random subset of those at i55 < 0, i.e. 
they have the same SFR distribution. The top panel of 
Fig. 3 shows that we cannot reject the null hypothesis: 
the sample with ^5 < 0 and the sample with ^5 > 0 are in¬ 
distinguishable. Similarly, a Pearson correlation test gives 
a P-value of 0.67. In other words, there could be a corre¬ 
lation only at the 32% conhdence level. We repeated the 
analysis for higher redshift slices, (zg = 3.1, 3.2, and 3.3) 
and Hnd similar results. 

Using the GIF2 simulations of section 3, we performed 
Monte-Garlo simulations to test whether or not the scatter 
in the Sg-pg relation is responsible for our null-result. We 
selected randomly about a third of the simulated LBGs 
to match our sample size of 334 LBGs. These simulated 
LBGs were assigned a SFR using SFRgim = (ps/ < ps > 

< SFR >, where p5 is the same as in Fig. 2, s is the 
slope of the simulated SFR-density relation, and < SFR > 
is the observed mean SFR in our data. Noise with the 
same properties as in the data (i.e. (t(SFR) ocSFR) was 
added to SFRsim. Note that these mock catalogs have the 
scatter between II5 and p5 (Fig. 2) built in. We repeated 
our analysis of SFRgim versus E5 1000 times for values of 
s spanning [—0.5,0.!]. 

At z = 0, the SFR decreases by a factor of 5 over 2 
decade in density (Gomez et al. 2003; Balogh et al. 2004), 
i.e. s ~ log(l/5)/2 ~ —0.35. In this case, the Pearson 
correlation test gives a probability smaller than 10“® (or 
~ 5 -(t) of Hnding no correlation. Even if we artihcially 
increase the noise by a factor as large as 2.5, as if we had 
underestimated our errors by such a factor, we Hnd that 
there is less than a 0.1% chance of finding no correlation if 
the z = 3 SFR-density relation were as steep as at z = 0. 
If the z = 3 SFR-density relation were not as steep, but 
were 3.6 times flatter (s = —0.07), we would still have 
detected a correlation at the 2 -a level. 

As advocated by Hogg et al. (2003), because the local 
density has a signal-to-noise ratio much lower than the 
other physical quantities, it may be preferable to compute 
the mean density at constant galaxy properties (e.g., SFR, 
M, color, etc.). The left panel of Fig. 4 shows the over¬ 
density Js as a function of SFR. The running median and 
the 1-CT spread are shown by the thick line and the error 
bars, respectively. Again, this plot reveals no evidence of 
an SFR-density relation. 

We also investigate the dependence of galaxy color on 
density. The right panel of Fig. 4 shows the over-density 
S 5 as a function of the observed-frame V — I. Again, the 
running median and the l-cr spread are shown. The rest- 
frame UV colors of LBGs do not appear to change with 
environment. We note that this non dependence of rest- 
frame UV color for blue galaxies such as LBGs is similar 
to the results of z = 0 surveys: in contrast to galaxies 
on the red sequence (e.g. Bell et al. 2004), the mean color 
of blue galaxies is only weakly dependent on environment 
(Blanton et al. 2003; Hogg et al. 2003, 2004; Balogh et al. 


2004). However, we cannot rule out an environment de¬ 
pendence of the mean rest-frame optical colors (or age) 
of LBGs given that we sample only Arest ^ 2000A. Fu¬ 
ture rest-frame optical observations should address possi¬ 
ble variations of the rest-frame optical colors of LBGs with 
environment. 

6. SUMMARY AND CONCLUSIONS 

From our wide-held images, covering a total of 0.90 
deg^, we selected z ~ 3 LBGs in several redshift slices 
100/1“^ Mpc (co-moving) deep spanning the redshift range 
2.9 < z < 3.4. We computed the SFR from the UV lu¬ 
minosity at rest wavelength 2000 A. Using mock catalogs 
from the GIF2 simulations, we show that our photomet¬ 
ric redshift accuracy is sufficient to statistically distinguish 
between low- and high-density regions using projected den¬ 
sity estimators such as the density within the fifth nearest 
neighbor S5. Our main results are as follows: we find 
that (1) there is no evidence of an SFR-density relation at 
z = 3, and (2) the rest-frame UV colors of LBGs do not 
appear to change with environment. 

Using Monte-Garlo simulations and the same projected 
density estimators that we applied to our data, we find 
a probability smaller than 10“® of finding no correlation 
if the steep z = 0 SFR-density relation were present in 
our data. If the z = 3 SFR-density relation were about 
four times flatter than the z = 0 relation, we would still 
have detected a correlation at the 2-a level. We conclude 
that, unless the z = 3 SFR-density relation is even flat¬ 
ter, we would have detected it in our data. This is in 
sharp contrast to surveys at z = 0 and z ~ 0.4 that have 
shown that the mean SFR is strongly dependent on the 
local galaxy density (e.g., Lewis et al. 2002; Gomez et al. 
2003; Balogh et al. 2004; Kodama et al. 2004). 

Very few predictions have been made of how the SFR- 
density relation should evolve from z ~ 3 to the present. 
Recently, Keres et al. (2004), using smooth particle hydro- 
dynamic cosmological simulations, showed that the SFR- 
density relation is expected to be present as early as z ~ 2, 
but not at z ~ 3 (see their Fig. 13). See Keres et al. (2004) 
and Birnboim & Dekel (2003) for a detailed description 
of the physical mechanisms at play. Our observed non¬ 
dependence of SFR on environment supports the theoret¬ 
ical results of Keres et al. (2004). 

Naturally, our sample is not representative of the en¬ 
tire galaxy population at z = 3; it is biased towards blue 
star-forming galaxies and does not include the red popu¬ 
lation unveiled by the Faint Infrared Extragalactic Survey 
(Labbe et al. 2003). With the present data, we are unable 
to rule out the existence of a SFR-density relation if one 
were to include this redder population. 
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Fig. 1.— The X—Y distribution (dots) of 50% (randomly selected) of the simulated galaxies of the GIF2 simulations (Gao et al. 2004). The 
filled triangles show all the simulated LBGs. Features of large-scale structure are visible in the simulated LBG population, considering that 
the galaxy positions are projected along the ^-direction over 100/i“^ Mpc. This is quantified in Fig 2. 
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Fig. 3.— Bottom: For LBGs selected to be at 2 phot = 3.0 ih 0.075 (see text), we plot the dust-corrected SFRq against the over-density ^5 
obtained from E 5 , the surface density within the fifth nearest neighbor. Error bars on SFR are shown for 20 (randomly selected) LBGs. Top: 
Results for the KS-test. For a given over-density ^ 5 , the probability that the distribution of galaxies with over-densities > ^5 is drawn from 
the reference distribution (defined as <55 < 0, left of the vertical dashed line) is shown. The distribution in the high density bins is statistically 
the same as in the low density bins (the 95% rejection level is indicated by the dotted line). 
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Fig. 4.— Left: The density-SFR relation at 2 = 3. The filled circles show the median over-density <^5 in bins of SFRq (bin size shown by the 
horizontal error bar). Right: Density-color relation at 2 = 3. The filled circles show the median over-density <^5 in color bins (bin size shown 
by the horizontal error bar). In both panels, the thick line shows the running median, and the vertical error bars the l-a spread, respectively. 
The colors and SFR of LBGs are independent of the density of the surrounding galaxy. 





